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Abstract 

The environment near super massive black holes (SMBHs) in galactic nuclei contains a large number 
of stars and compact objects. A fraction of these are likely to be members of binaries. Here we discuss 
the binary population of stellar black holes and neutron stars near SMBHs and focus on the secular 
evolution of such binaries, due to the perturbation by the SMBH. Binaries with highly inclined orbits 
in respect to their orbit around the SMBH are strongly affected by secular Kozai processes, which 
periodically change their eccentricities and inclinations (Kozai-cycles). During periapsis approach, 
at the highest eccentricities during the Kozai-cycles, gravitational wave emission becomes highly 
efficient. Some binaries in this environment can inspiral and coalesce at timescales much shorter 
than a Hubble time and much shorter than similar binaries which do not reside near a SMBH. The 
close environment of SMBHs could therefore serve as catalyst for the inspiral and coalescence of 
binaries, and strongly affect their orbital properties. Such compact binaries would be detectable 
as gravitational wave (GW) sources by the next generation of GW detectors (e.g. advanced- LIGO). 
About 0.5 % of such nuclear merging binaries will enter the LIGO observational window while on orbit 
that are still very eccentric (e > 0.5). The efficient gravitational wave analysis for such systems would 
therefore require the use of eccentric templates. We also find that binaries very close to the SMBH 
could evolve through a complex dynamical (non-secular) evolution leading to emission of several GW 
pulses during only a few yrs (though these are likely to be rare). Finally, we note that the formation 
of close stellar binaries, X-ray binaries and their merger products could be induced by similar secular 
processes, combined with tidal friction rather than GW emission as in the case of compact object 
binaries. 
Subject headings: Gravitational waves - binaries - Galaxies: Milky Way galaxy. 



I. INTRODUCTION 

The majority of stars reside in binary and higher mul- 
tiplic ity systems (jRaghavan et al.l l2010t iRemage Evans! 
1201 It ) , and are observed to exist both in the field as well 
as in dense stellar environments such as globular clusters 
and galactic nuclei. In dense clusters the dynamical evo- 
lution of binaries can be strongly affected by encounters 
with others stars. In the field, encounters with other 
stars are rare, and do not affect the binary evolution. 
Such binaries may still dynamically evolve if they re- 
side in triple (or higher multiplicity) systems. The inner 
binary in a triple can be affected by the perturbation 
of the third, outer companion, and evolve through long 
term secular processes. In particular, triples in which the 
relative inclination between the orbit of the inner binary 
and the outer bin ary is high ca n evolve th rough Kozai- 
Lidov oscillations (jKoza i 1962; L idovl[l962n in which the 
eccentricity and inclination of the inner binary signifi- 
cantly change in a periodic/quasi-periodic manner. 

In this study we explore the evolution of binaries in 
the environments of SMBHs which are thought to re- 
side at the center of most galaxies (e.g. iGebhardt et al.l 
120001 12003| ) ; throughout this paper we focus on galactic 
nuclei similar to the Galactic center (GC) of our own 
galaxy. Close to the center, where the gravitational po- 
tential is dominated by the mass of the SMBH, binaries 
are bound to the SMBH and effectively form a triple 



system, in which the binary orbit around the SMBH is 
the outer orbit of the triple. Binaries near a SMBH may 
therefore evolve both through gravitational scattering by 
other stars in the dense nuclear cluster, as well as through 
secul ar evolution in the SMBH-stellar binary triple sys- 
tem (jAntonini et al.l[201Cl) . Various types of binaries ex- 
ist in such environments, including both main-sequence 
(MS) binaries, compact object (CO) such as white dwarf 
(WDs), neutron star (NS) or stellar black hole (BH) bi- 
naries or mixed CO-MS binaries. In this paper we focus 
on CO-binaries, their evolution into short periods and 
the way in which they may become potentially observ- 
able GW sources for the advanced laser interferometer 
GW observatory (LIGO). Scattering with other stars is 
also accounted for, but only in the sense of the lifetime 
of binaries in the GC, i.e. we do not follow the evolution 
of binaries which are disrupted (evaporated) through en- 
counters with other stars at timescales shorter than the 
relevant secular timescales (Kozai evolution). 

In the following we begin with a discussion of bina- 
ries in the GC and their expected properties (Section 
[2]), we then discuss the relevant timescales for the evolu- 
tion of stars and binaries in galactic nuclei (Section [S]). 
In Section |4] we discuss the secular evolution of binaries 
near a SMBH, and their potential coalescence due to a 
combined Kozai cycles-GW loss mechanism; we then cal- 
culate the type and properties of the GW sources which 
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form through this process. The secular evohition of other 
types of binaries is briefly discussed in Section [5] . We 
summarize our results in Section [H] 

2. BINARIES IN THE GALACTIC CENTER 

At present, little is known about binaries in the 
GC. Massive binaries can be currentl y detected only 
by observati ons of stellar variabil ity lO tt ct alJ 119991: 
DePov et alJ l200l iMartins et all mm Rafelski et all 



20071 1 ■ which can mainly detect eclips ing; or ellipsoidal 
variable binaries (jRafelski et alJ 120071 ). The fraction 
of such bi naries among 0-star s outside the GC i s 
2 - 11% fGarmanveiaP flOSOt IRafelski et aP l2007l) . 
jRafelski et al. (2007|) conducted a variability survey of 
the central 5" x 5" o f the GC and found one such bi- 
nary (also observed bvlOtt et al.lll999l:lDePov et al.l l2004l: 
IMartins et alll2006l : IRafelski et al.ll2007| ) among 15 mas- 
sive stars (7%); a fraction consistent with the Galac- 
tic one. Other binar y candidates were a lso found: a 
colliding wind bi narv (IRafelski et al.|[2007ll and another 
eclipsing binary (jPeeples et al.l 120071 ). Although these 
data are still insufficient for drawing strong conclusions 
about the total massive binary fraction in the GC, they 
do suggest that it is comparable to that observed out- 
side the GC, among massive binaries. For reference, 
the massive binary f raction in t he Solar neighb orhood 
is very high (> 70%:lGarmanv e t al. 1980; Maso n et all 
199a iKouwenhoven et all l2007t .Remage Evansi 120111 : 



Kobulnicky fc Fryeil 120071 ) , and most massive binaries 
have semi-major axes of up to a few AU. Given our 
poor understanding of binary formation, it is unknown 
whether star formation in the unusual GC environment 
results in similar binary properties. 

No observational data exist regarding low mass bina- 
ries in the GC. Little is known about compact binaries 
in the GC, though the obs erved overabundan ce of X- 
ray sources in the central pc (jMuno et al.ll2005l ) suggests 
they are not rare but their properties and origin are not 
known. 

We focus on binaries very close to the SMBH, 
where secular Kozai evolution affects the binary dynam- 
ics. However, this population is continuously destroyed 
through various processes (evaporation, mergers, disrup- 
tion). We therefore first discuss the binary population 
outside the central region near the SMBH (>~ 2pc), 
which serves as a continuous source term repopulating 
the binaries in the central region. We then discuss the 
processes which can repopulate binaries in the central 
region, both from populations outside and from in-situ 
formation of new binaries. 

2.1. Binaries outside the radius of influence of the 
SMBH 

In the following we discuss the binary population out- 
side the central region near the SMBH (>^ 2pc). This 
binary population serves as a source term for binaries in 
the central region on which we focus, which are continu- 
ously destroyed through various processes. 

The binary fraction and semi-major axis distribution 
of observed main sequence stars in the solar neighbor- 
hood have been analyzed i n many binar y surveys (e.g. 
iDuauennov fc Mavoil (|1991l ): lLadal(|2006D and references 
within). For compact objects the observational data is 



still lacking. However, many theoretical analysis have 
been done using stellar evolution and binary synthesis 
codes. Most of these studies focused on close binaries, 
usually in the context of CVs, x-ray binaries or gravita- 
tional waves sources (which are relatively short lived and 
compose only a very small fraction of the binary pop- 
ulation), and did not analyze the general distribution, 
which is needed for our analysis. Nevertheless, the gen- 
eral semi-major axis distribution of COs binaries have 
been found for isolated (i.e. not taking into account 
dynamical effects which could be im portant in dense 
stella r environments) WD- MS binaries (iWil ems fc Kolbj 
I2004D and BH binaries (Belczvnski et al.l 12004^ Ebr 
NSs, one can use the observed pulsar binaries 
population as fo und in t he ATNF pulsar catalog 
(at http : / / www, atnf . cs iro. au /research /pulsar /psrcat / ; 
iManchester et all 120051 ). The initial period distribution 
we consider for the different binary populations is shown 
in Fig. [T] These distributions do not take into account 
dynamical evolution of binaries in the dense environment 
near the SMBH, mainly impor tant in evapo rating wide 
binaries (so called soft iDinaries; lHeggiellT975[ ) due to dis- 
tant stellar encounters. They also do not account for 
the long term evolution of close compact object binaries 
due to GW emission. We account for these evolutionary 
aspects using simplified methods. 

The primordial binary fraction (not taking into ac- 
count dynamical evolution such as evaporation of the 
binaries due to encounters), fbm , is different for the dif- 
ferent binary types, MS, WD, NS or BH binaries. It de- 
pends both on the progenitors binary fraction, /^™ (i.e. 
BHs and NSs progenitors are high mass stellar binaries 
whereas WD binary progenitors have lower masses) and 
fraction of surviving binaries, fsurv, after going through 
their evolutionary stages . 

The binary fraction of massive stars (/^™ — 0.7) is 
more than two times that of lower mass stars (/^™ = 0.3; 
iKobulnickv fc Frveil 12001 iLada 2006). Stellar evolution 
further changes these fractions. Thus, for obtaining /(,„ 
we either take the observed binary fractions in the So- 
lar neighborhood (in the cases of MS binaries and NS 
binaries) or use the theoretical models for the surviv- 
ability of the different bina ries (jWillems fc Kolbl I2002L 
12004 iBelczynski et aDl2004[ ) to find fsurv, and multiply 
it by the progenitors binary fractions f^[° (i.e. fbin = 
funfsurv). Wc finally get fun = 0.7, 0.3, 0.3, 0.07, 0.1 
(for high and low mass MS stars, WDs, NSs and BHs 
respectively) as the evolution of lower mass stars (and 
WDs) changes the semi-major axis distribution of the 
binaries but does not result in their disruption, whereas 
the evolution of BHs and NSs results in most cases in the 
binary disruption. We adopt these initial binary frac- 
tions throughout our calculations. The fractions of CO 
binaries used here already account for the initial high 
binary fraction and the stellar evolution of these bina- 
ries. However, these binary fractions are also expected 
to be reduced due to dynamical encounters with other 
stars and through coalescence due to GW emission (for 
compact objects) and tidal friction. In the following cal- 
culations (specifically in section H75|) we assume the ini- 
tial CO binary fractions {fbin) in the GC are similar to 
those obtained here, where the effects due to dynami- 
cal encounters and Kozai evolution of compact binaries 
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Fig. 1. — The period distribution of primordial bina- 
ries (without dynamical effects). Main sequence stars 
(solid thick line) di s tribut ion taken according to observations 
IDuquennov fc Mayor! 1 119911 ). WD distribution (thin dashed line) 
is t aken according to the theoretical model of IWillems & Kolbl 
I|20d4 ). NS binaries distribution (thin solid line) is taken 
according to distribution of binary pulsars in the ATN F 
database (a t http://www.atnf.csiro.au/research/pulsar/ps rcat/| ; 
IMa nchcster et al. ^UUEj). BH binaries distribution (thick triple 
peaked dashed line) is ta ken according to the theoretical model 
of rSelczynski et al.. (,2001) . 

are accounted for in a simplified manner, as we discussed 
later on. Evolution due to tidal frction which is impor- 
tant for MS stars will be discussed elsewhere, and is only 
shortly discussed here. 

2.2. The binary fraction near a SMBH 

As discussed above, the binary fraction in the GC is 
unknown; and the current observational data is insuffi- 
cient to prove any meaningful constraints. Very few the- 
oretical studies explored the role and e volution of bina- 
ries near SMBH. lBahcah fc WoJJ ()1977[ ) briefly discussed 
the effect of binaries on the dynamics of stars in a cusp 
around a SMBH. They showed that their effect on the 
cusp dynamics is small, since the vast majority of these 
binaries are soft, and the total energy which they can 
transfer to the cusp stars is small in respec t to their or- 
bital e nergy (in the potential of the SMBHl. lPerets et al.l 
(j2008| ) suggested that binaries in the observed stellar disk 
in the GC can have a non-ne gligible effect on th e evolu- 
tion of the cold disk (see also lCuadra et al]l2008[ ). Bina- 
ries can also be important for a fast deposit ion of stars 
close to a SMBH through binary disruption (iHillsl 119881: 



Gould fcQuillenll200a iPerets et al.|[2UOT: iMadigan et aLl 
2009( 1. These studies did not discuss the binary popula- 
tion very close to the SMBH; the only study to explore 
the general popu lation of binaries near SMBH was done 
by iHopmanI (|2069l) . even though this pioneering study 
accounted only for regular relaxation processes. A full 
study of the binary population near SMBH is highly de- 
sirable. However, this is beyond the scope of this study 
which mainly focuses on the secular evolution of binaries 
due to the perturbation by a SMBH. Here we assume a 
simple model in which the binary population in galac- 
tic nuclei follows their population in the field, but we 
also account for their evaporation through taking their 
effective lifetimes to be the evaporation timescale at the 



distance from the SMBH where they reside. We discuss 
several possible replenishment of binaries in Section 14.41 

3. TIME SCALES 

In the dense stellar environment of a galactic nucleus, 
otherwise rare dynamical processes can take place and 
affect the evolution of the stellar population near the 
center. Here we compute the time-scale associated with 
such processes in the Galactic center and compare it to 
the Kozai time-scale associated with the secular evolu- 
tion of binaries due to the perturbations by a central 
SMBH. 

We consider stellar-mass binaries with components 
of masses mg and mi orbiting a super-massive black 
hole (SMBH) of mass M,. We denote the eccentricities 
of the inner and outer binary, respectively, as ei and 62, 
and semi-major axes ai and 02. We also define gi as the 
argument of periapsis of the inner binary relative to the 
line of the descending node, and / as the mutual orbital 
inclination of the inner binary with respect to the outer 
orbit. 

The mass density near the SMBH was modeled using 
the broken power-law density profile: 



P{r) = Po 



ro 



1 



(7-1.8)/2 



(1) 



with 7 the inner slope index and tq = 0.5pc. Setting 
Po = 5.2 X IO^Mqpc"'^ this gives a good fit to the ob- 
served space density at t he Galactic center o utside the 
core, normalized at 1 pc (jSchodel et al.|[2009D . For the 
slope of the inner density profile we assumed two differ- 
ent values: 7 = 0.5, representative of t he observed dis- 
tribution of stars at the Galactic center ( Buchholz et al.l 
[2009; Do et al. 2009; Bartko et al. 2010), and 7 = 1.8 
corresponding approximately to a nearly-relaxed config- 
uration of stars n ear a dominating point-mass potential 
(|Alexandeill20050 . We set M, = 4 x 10*^ Mm similar to 
the mass o f the SMBH at th e Galactic center (jGhez et al.l 
2008; Gillc ssen et aTI[200l . 

Kozai timescale: The Kozai time-scale can be written 
in terms of the masses of the three bodies, the eccentric- 
ity of the outer binary and their semi-major axes as: 
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where Pi is the period of the inner binary, P2 the period 
of the outer orbit, and Mh = mo + mi. 

Extra sources of periapsis precession: Kozai Cycles can 
be suppressed by additional sources of apsidal precession. 
The most relevant processes to consider are relativistic 
precession and, in the case of stellar binaries, precession 
due to the tidal bulge raised on each of the inner bi- 
nary's components. Writing the timescale of relativistic 
precession as: 
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Fig. 2. — Time scales at the Galactic center (dashed lines) compared to the Kozai period (Tk; continue lines) for a (10 + 10) A/q (red 
lines) and (1 + 1) Mq (black lines) binaries with ai = 1 AU. Rectangular regions show the age and radial extend of the S-stars and disk 
of WR/O stars. Dotted lines give the radii at which Kozai cycles are fully suppressed by relativistic precession in the inner binary. Open 
circles indicate the radius where the time-scale of periapsis precession due to tidal bulges equals Tk assuming an inner eccentricity e\ such 
that the closest approach between the stars is twice the sum of their radii, i.e., ei = 1 — 2 X (i?o + ^i) jo-l- 



with G the gravitational constant and c the speed of 
hght, it can be shown that general relativistic preces- 
sion in the inne r binary stops the Kozai oscillations when 
(jHollvwood fc Meha 1997 : Blaes et al. 2002) : 
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Tidal bulges always tend to promote periapsis pre- 
cession and therefore suppress Kozai cycles. The time 
sca le of periapsis precess i on due to non dissipative tides 
is (IKiseleva et all 119981 : lEggleton fc Kiseleva-EggletonI 
[200llh 
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with Rq and i?i the stellar radii and ko and fci the tidal 
Love numbers. Assuming equal objects in the binary and 
ko = ki = 0.01, we find: 



T 



TidG = 



1.3 X 10^ V 



/ Ql 

VlA1 



lAU 



13/2 / M, 



2Mr. 



-1/2 



i?0 
IRr. 



{l-ei) 



2\5 



8 + 12ef -I- ef 



-5 

(6) 
(7) 



where the stellar rad ii are obtained using the expression 
(jHansenet al.ll2004D : 



R — Rr. 



M 

Mo 



0.75 



(8) 



Comparing Equations ([3]) and ^ it can be seen that in 
general GR precession is initially the most rapid mech- 
anism inducing apsidal motion, but, for very large ec- 
centricities, when, for instance, the system approaches 
the maximum of a Kozai cycle, one expects tides to be- 
come important to the evolution. Consequently, for stel- 
lar binaries with relatively large semi-major axes, very 



large eccentricities (eventually resulting in mass trans- 
fer or even collisions) could be avoided due to apsidal 
precession induced by tidal bulges. 

Binary evaporation: In a dense environment, binaries 
may evaporate due to dynamical interactions with field 
stars if 

\E\/{M,a^) < 1, (9) 

with E the internal orbital energy of the binary, and 
a the one-dimensional velocity dispersion of the stel- 
lar background. The evap oration timescale is given by 
(jBinnev fc Tremaind [19871) : 
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where InA is the Coulomb logarithm, a and p are the ID 
velocity dispersion and mass-density, and to is the mass 
of a typical object in the system assumed to be solar 
mass stars in what follows. The velocity dispersion, a, is 
calculated from the Jeans equation: 

/•OO 

p{r)cr{rf ^G dr'r'-^ [M, + M^{< r')] p{r'), (11) 



with Mi,(< r) the total mass in stars within r. 

Two-body relaxation: Assuming equal-mass stars and 
an isotropic velocity distribution, the local t wo- body 
(non resonant) relaxation time is defined as (jSpitzed 
[19871) : 



TNR^l.ex 10i"yr 
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306 km s-i 
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m , , , 



2.1 X IO6M0PC-3 

Setting 7 — 0.5 in the density model of Equation (|T]), this 
time results somewhat larger than the age of the Galaxy 
[r^ 10^'^ yr) when com puted at the radius of influence 
(- 3 pc) of the SMBH (|Merritti[20Tdh . 

Resonant relaxation: In the dense stellar environment 
near a SMBH, as long as the relativistic precession time 



Evolution of binaries near massive black holes 



scale is much longer than the orbital period, the mech- 
anism that dominates the scattering of stars onto high- 
eccentricity orbits is resonant relaxation. Because in the 
potential of a point-mass the orbits are fixed ellipses, 
perturbations on a test particle are not random but cor- 
related. The residual torque |Q| sa ^/NGm/r, exerted by 
the N randomly oriented, orbit-averaged mass distribu- 
tions of the surrounding stars, induces coherent changes 
in angular momentum AL = Qt on timescales t < Tcoh- 
Where the coherent time tcoh is fixed by the mechanism 
that most rapidly causes the orbits to precess (e.g, mass 
precession, relativistic precession). The angular momen- 
tum relaxation time associated with resonant relaxation 
is: 



where Lc = y/GM^a is the angular momentum of an 
orbit with radius r '^ a and jALco/il ^ |QT'co/j| is the 
accumulated change over the coherence time. Assum- 
ing that the precession is determined by the mean field 
of stars (i.e., mass precession) , the angular momentum 
relaxation time is (iRauch fc Tremainelll996l ): 
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(13) 
At small radii relativistic precession becomes efficient 
at suppressing resonant relaxation and the gravitational 
perturbations are dominated by classical (two-body) 
non-resonant relaxation. The value of the angular mo- 
mentum (the S chwarzschild barrier ) at which this tran- 
sition occurs is (jMerritt et al.l[2011l ): 
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where N is the number of stars within radius a2 and Csb 
is a constant of order of unity. 

Vector resonant relaxation: An additional fundamen- 
tal process in changing the external orbit of the binary 
is vector resonant relaxation (VRR). The external or- 
bital plane is randomized non-c oherently on the vec- 
tor res onant relaxation timescale (jHopman fc Alexandeii 
l2006aD : 
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VRR changes the direction but not the magnitude of 
the external orbit angular momentum. The randomiza- 
tion process of the external orbital plane can possibly 
produce variations in the mutual inclination of the bi- 
nary orbit with respect to its orbit around the SMBH. 

In Figure 2 these different timescales are compared for 
the case of a circular orbit and for an eccentric orbit 
with 62 = 0.95. In this latter case the corresponding 
timescales are obtained by multiplying the quantities of 



Equations (ini),(II31) and ^5^ by a factor (1-62). To com- 
pensate for the fact that when the external eccentricity 
is large the binary spends most of the time at apoapsis, 
we divided the evaporation time given in Equation (|10|) 



by a factor ^1 — e|. At the radii of interest, for either 
the cusp and the core model, the Kozai timescale re- 
sults much shorter than any other timescale associated 
with perturbations from dynamical interaction with field 
stars. VRR is the most rapid mechanism associated with 
gravitational interaction with field stars. For 02 < 0.1 pc, 
since Trr^v ^ Tk, VRR might lead to an increase of the 
number of binaries undergoing Kozai resonances with the 
SMBH: large inclined binaries will have already experi- 
enced multiple Kozai cycles and eventually merge before 
VRR could change appreciably their orbits, while low in- 
clined binaries could move to high inclinations via VRR 
and enter the Kozai regime. The degree to which the 
binary dynamics depend on VRR and other dynamical 
processes is a complicated topic that would most likely 
require direct A'^-body simulations describing dense clus- 
ters of stars around SMBHs. In what follows, as a first 
approximation, we neglect VRR and use the time for bi- 
nary evaporation (Tev) as the timescale of reference to 
define the parameter space where our integrations are 
valid. 

From Figure 2 we see that within ^ 0.5 pc, i.e., within 
the limiting radius of the stellar disk at the Galactic cen- 
ter, most of the binaries with initially high inclined orbits 
would have performed multiple Kozai cycles during their 
life time. At ^ 0.2 pc, Kozai resonances are completely 
suppressed by relativistic precession in the inner binary. 
However, since increasing the external orbital eccentric- 
ity reduces the Kozai period but leaves the timescale for 
relativistic precession unchanged, for a given orbital ra- 
dius it is always possible to choose 62 in order to have 
Tk » Tgr- For instance, setting 62 > 0.9, relativistic 
precession becomes important only at radii of order par- 
sec (middle panel in the figure). However, very large ec- 
centricities could result in a very short distance of closest 
approach of the binary to the SMBH, possibly resulting 
in the tidal break up of the system. This occurs when 
the binary approa ches the SMBH within at a distance 
(jMilleret al.ll2005D : 
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In this section, we have assumed that the perturbers 
are solar mass objects, but, we note here that gravi- 
tational scattering can be dominated by stellar black 
holes. This can be the case if the age of the nucleus 
is larger than the mass-segregation time scale of the 
heavy component. Assuming standard initial mass func- 
tions, in a relaxed model the mass i n black holes within 
0.01 PC from SgrA* is A^rh ?^ 10^ (JFreitag et al.l 12001 
iHopman fc Alexandedl2006ah and Gravitational pertur- 
bations will be dominated by stellar black holes. In this 
case the time scale for evaporation will be somewhat 
shorten than that implied in Figure 2 within 1 mpc, but it 
will be essentially unaffected at larger radii. In the case of 
a shallow density cusp around a massive black hole (i.e.. 
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Fig. 3. — Time evolution of the gravitational wave frequency (/gw)i eccentricity (ei) and semi-major axis of binary black holes orbiting 
a SMBH. In these integrations the distance of closest approach of the binary to the SMBH is rper = 3.3 X r^t, with r^t the tidal disruption 
radius of the binary. We set the initial parameters to ei = 0, mo = mi = 10 Mq, I = 88°, gi = 0° and adopt different values for the orbital 
eccentricity e2 and inner semi-major axis ai. The nominal merger timescales of the binaries in isolation are ~ 10^'* yr (ai = 1 AU; bottom 
right panel), ~ 10^'' yr (ai = 0.5 AU; two bottom left most panels), ~ 10^" yr (ai = 0.1 AU; middle panels) ~ 10^ yr (ai = 0.05 AU; top 
panels). The presence of an external perturber, the SMBH in this case, causes periodic variation of the inner binary eccentricity greatly 
reducing the merger time-scales of these systems. 



7 < 0.5) dynamical friction falls essentially to zero inside 
the core due to the absence of stars moving more slowly 
than the local circular velocity (jAntonini fc Merritt I 
|2012| ) . This results in much lower densities of stellar rem- 
nants at small radii with respect to c ollisionally relaxed 
models (|Gualandris fc Merritt I I2012D . The distributed 
mass in this case will be therefore dominated by stars 
at all radii, and gravitational scattering from black holes 
can be neglected. 

4. GRAVITATIONAL WAVE SOURCES 

In the following we address the possibility that the 
merger time of compact binaries residing in galactic nu- 



clei can be significantly accelerated due to their interac- 
tion with a central SMBH. Some systems that otherwise 
would have no chance of merging before being disrupted 
by gravitational encounters with background stars (or 
stellar remnants) will actually merge due to combined 
effect of Kozai cycles and gravitational wave loss. For an 
isolated binary, the gravi tational wav e merger timescale 
can be approximated as (|Peterslll964l) : 
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Fig. 4. — Evolution onto the ei — /gw and ai — ei planes of 
a stellar black hole binary orbiting a SMBH on a circular orbit 
of radius 20 AU with / = 88°. This system enters the LIGO 
observational band with e ~ 1 and merges in ~ 20 yr after many 
periodic variations (~ 100) of the inner eccentricity have already 
occurred. 

Given the strong dependence of To on the eccentricity, 
the induced oscillations in the binary orbital eccentric- 
ity via secular Kozai processes can cle arly decrease dra - 
matically the merger timescale (e.g. lThonipsonll2Qll[ ). 
The GW signal emitted by such compact binaries will be 
completely dominated by repeated pulses emitted dur- 
ing periapsis passages and at hi gh eccentrici ties where 
the GW frequency is maximized (jGouldll201lD . 

4.1. N-body integrations 

In this section we describe high accuracy A^-body simu- 
lations that were used to study the evolution of a limited 
selection of systems representing binaries of BHs orbiting 
a SMBH. These simulations were carried out using the 
AR-CHAIN integrator (Mikkola & Mcrritt 2008), which 
includes PN corrections to all pairwise forces up to order 
PN2.5. The code employs an algorithmically regular- 
ized chain structure and the time-transformed leapfrog 
scheme to accurately trace the motion of tight binaries 
with arbitrarily large mass ratios. We considered circu- 
lar and equal mass binaries with components of masses 
10 Mq. The periapsis distance of the external binary 
orbit was 3.3 x rbt, and we fixed the initial inclination of 
the inner orbits with respect to the external orbit around 
the SMBH at / = 88°. We then varied 62 and oi in or- 
der to study the dependence of the merger timescale on 
these latter quantities. The choice of a large inclination 
was also motivated by the short merger timescale of such 
systems allowing for an efficient computation of their dy- 
namical evolution. 

Figure [3] gives the time evolution of semi-major axis, 
eccentricity and gravitational wave frequency (/gw) of 



the inner binaries. The pea k gravitat ional wave fre- 
quency was approximated lay (jWenl [20031 ): 
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with semi-major axis and eccentricity obtained from the 
relative r adial expressions given in E quations (3.6 a) and 
(3.6 b) of IDamour fc Deruelld (IT985I) . 

The oscillations in the binary inner eccentricity and 
gravitational wave frequency induced by the SMBH per- 
turbations speed up the merger time by many orders of 
magnitude. For example, the nominal gravitational wave 
merger timescale for the cases shown in the middle pan- 
els of the figure without SMBH is ~ 10^° yr, while it is 
just few years in our simulations. The oscillations in 
eccentricity and inclination observed in the computed 
orbits are reminiscent of the "Kozai cycles" generally 
discussed in the context of hierarchi cal triple systems 
(e.g.lHarringtonll968l:lWu et al. l20"07HNaoz et al.i2011al: 
iProdan fc Murravll2011[ ). However, the distance of clos- 
est approach to the SMBH for the orbits of Figure [3] 
is too short (about 3 times the binary tidal disruption 
radius) for the system to be treated exactly as hierar- 
chical. In other words, the strong interaction at periap- 
sis causes high-order terms (beyond the octupole-order) 
in the equations of motion to become important to the 
evolution. The amplitude of the eccentricity oscillations 
changes essentially at each cycle, eventually making the 
inner binary periapsis distance small enough that effi- 
cient GW loss takes place. Interesting, we found that the 
changes in eccentricity, for initially high eccentric orbits 
(^2 ^ 0.5), always occur very rapidly during the closest 
approach of the binary to the SMBH rather than in a 
smooth-like way along the entire orbit as it would be the 
case if we had used the secular perturbation equations in- 
stead of direct A'^-body integrations. This behavior is ev- 
ident in Figure 131 but see also Figure 7 in lAntonini et al.l 

dMa). 

We note also that in some cases (e.g., upper left panel 
of Figure [3|) it is possible that the binary first enters a 
regime in which GW starts to dominate the evolution, 
suddenly causing a drop of the binary semi-major axis; 
during this phase the binary is then brought back to a 
small-eccentricity configuration where GW loss becomes 
again negligible with the subsequent evolution occurring 
at approximately constant semi-major axis. 

Figure[4]displays the evolution of a stellar black hole bi- 
nary with ai =0.1 AU, inclination / = 88° and gi = 0°. 
The binary is on a circular orbit of radius 20 AU from 
the central black hole. During the oscillations in eccen- 
tricities, at the maximum of a Kozai cycle, the binary 
gravitational wave frequency is such that system could 
be a detectable source for LISA-like gravitational wave 
detectors. The signal from such sources in a LISA-like 
detector would result in repeated pulses emitted during 
periapsis approach. The system merges in only ~ 20 yr, 
after many eccentricity oscillations have occurred and en- 
ters the LIGO observational window (> lOHz) with a 
very large eccentricity (ei « 0.99995) before the fast or - 
bital circularization due to GW loss begins. iWenI (|2003[ ) 
showed that, in the contest of stellar black hole triples 
that might form in the core of globular clusters, for some 
particular configurations it is possible that the inner bi- 



Antonini & Perets 




1000 2000 3000 4000 



1000 2000 3000 4000 



1000 2000 3000 4000 



2/a 



Fig. 5. — Merger time (continue lines) and evaporation time (dot-dashed lines) of stellar black hole binaries near a SMBH. The evaporation 
time was computed under the assumption of either a cusp or a core in the density of stars near the center. Highly inclined systems close to 
the central black hole merge in a time typically much shorter than the timescale over which they would evaporate due to close encounters 
with field stars. Dashed lines give the inspiral time of the compact binary into the SMBH due to GW loss. Here we used mg = mi = 10 Mq, 
ei = 62 = 0.1, and gi =0°. Black filled circles were obtained by direct 3-body integrations including PN corrections up to order 2.5. 



nary GW frequency enters the LIGO observational band 
before the maximum eccentricity expected for the Kozai 
oscillations is reached. In these cases the eccentricities 
can remain significantly large {^ 0.9) at /gw ^ lOHz. 
We stress that the system shown in Figure |4] is quali- 
tatively different from those discussed in this previous 
work: only after multiple periodic oscillations in the or- 
bital elements have already occurred, the binary reaches 
very high eccentricities and enter the LIGO frequency 
band. Such a situation occurs before gravitational radi- 
ation becomes important, allowing the eccentricity to be 
so large at such high frequencies. 

4.2. Approximate method 

Given the significant computational resources per run 
made necessary by the high accuracy of the iV-body in- 
tegrations presented in the previous section, no attempt 
was made to systematically explore the large parameter 
space of the three body interactions. In order to per- 
form such broader exploration we use in what follows an 
approximate method which is detailed below. 

If the external orbit of the inner binary is wide enough, 
i.e., the gravitational perturbations from the SMBH on 
the binary are weak, the system can be regarded as a 
hierarchical triple and it is possible to treat the dynamics 
of the entire system as if it consists of an inner binary 
of point masses uiq and mi and an external binary of 
masses Mi, = mo + mi and Af,. We define the angular 
momenta Gi and G2 of the inner and outer binary and 
total angular momentum H = Gi + G2: 



and 



Gi = TOoTOi 



G2 = Mb M, 
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Mb + M, 
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The quadrupole-level secular perturbation equations 
describing the evolution of the inner binary orbital el- 



ements are (|Ford et al.ll2000[ ) : 
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and 
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Equations (22)- (25) can be considered as describing the 
interaction between two weighted ellipses instead of point 
masses in orbits and are accurate as long as the energy 
of the system and its angular momentum are approx- 
imately conserved quantities within each Kozai-cycle. 
These equations include the effects of relativistic pre- 
cession and energy loss due to gravitational wave emis- 
sion and they can be integrated numerically much more 



Evolution of binaries near massive black holes 



q 
d 



■To=1oV 25 AU 
a, =0.05 AU / 


10 AU 


: To=1.5x1o'°yr/ 
; a, =0.1 AU / 

/so d 


' i'„»r=5 AU - 


1 ' 1 
/ ' / 


20 Au/ ; 

/ / 
/ / 


: / ' / 

- ' 1 
/ / / 


/' / : 


^~^^~y'^^—^~,£C^ - -' 


1 r 


1 



l2011bH : 



10 100 1000 

(1+62)7(1-62) 




1-e, 

Fig. 6. — Upper panel: merger time (in units of the merger time of 
the binary evolved in isolation) versus rapo/''per = (l+e2)/(l — 62)- 
At large semi- major axis, relativistic precession "de-tunes" the sec- 
ular Kozai effects and decreases the maximum eccentricity achieved 
during a cycle. This results in longer merger times. Bottom panel: 
merger time as a function of the internal eccentricity of the com- 
pact binary. For e\ > 0.7, the effect of the Kozai resonance on the 
merge time becomes progressively less important as e\ increases. 
In these computations we used mo = m\ = 10 Mq, ei = 0.1, 
/ = 85° and gi =0°. 



rapidly than the full equations of motion. They are used 
here to investigate the dynamical evolution of compact 
binaries near massive black holes. 

We set the mass of the binary components similar 
to that typical of stellar black holes: mo = rrii — 
10 Mo(Woosiey_e_t al. 2002). In this case Equations (22)- 
(24) are equivalent to the octupole-lev e l secu lar pertur- 
bation equations given in iBlaes et al.l (120021 ) . In fact, 
from Equation (19) of IBlaes et al.l (|2002f ) we see that 
when TTT-o = mi the octuple terms vanish. We also run 
integrations for unequal mass binaries and including oc- 
tuple terms to see the effect that higher order terms have 
on the binary evolution. These integrations gave results 
very similar to those obtained using the quadruple ap- 
proximation. The fact that the octuple terms have only 
a small effect on the binary dynamics is easily explained 
if we parametrize the relative contribution of th e octu- 
ple level approximation using the parameter (Naoz et al.l 
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(27) 



In our systems a2 ^ ai and even for very large eccen- 
tricities (e2 ~ 0.9) the contribution of the octuple terms 
is negligible. 

In Figure[5]we give the merger time, Tj/, as a function 
of the semi-major axis of the outer orbit, this last quan- 
tity is given in units of the initial binary semi-major axis 
(ai). We set ei = 62 = 0.1, ai = (0.1, 0.05, 0.01) AU 
and we consider several values of the inclination angle. 
The merger time is compared to the evaporation time 
(Tev) of the binary systems which we computed us- 
ing Equation ()10p . The time for evaporation in a shal- 
low density profile is greatly increased with respect to 
the cusp model and it is independent on radius since, 
close to the SMBH, both p and cr in Equation (|10l) go 
like ^ r^'^-^ . 

Highly inclined binaries (/ > 70°) that lie close to the 
central black hole inspiral and coalesce over timescales 
much shorter than a Hubble time and typically shorter 
than their evaporation timescale on such environments. 
At large distances from the SMBH relativistic preces- 
sion becomes efficient at suppressing Kozai cycles and 
the merger time approaches its value when the binary 
evolves in isolation. Notice that, since the perturber in 
our simulations is much more massive than the inner bi- 
nary, there is essentially no difference between retrograde 
and prograde orbits with the minimum time for merger 
occurring at / = 90°. Filled black circles in the figure 
were obtained by direct integrations of the three body 
system carried out using the ARCHAIN integrator. The 
agreement between direct integrations and the results 
obtained by using the secular perturbation equations is 
indeed very good. 

Dashed lines in Figure [5] give the time required for 
the stellar black holes to spiral into the central black 
hole. This timescale was approximated by considering 
the compact binary as a single point mass and by using 
iPeterd (|1964[ ) formula: 



T^,' 



:1.8 X 10^ V 



/20AfQ\ /4x IO^Mq 



V Mb 



M, 



(28) 



^2 

mpc 



)\i-eiy^' 



For small inclinations (< 70°), within 02/01 < 700, the 
merger time of the binary into the SMBH can be shorter 
than the time required for the stellar black holes to merge 
due to gravitational wave loss. However, the compact bi- 
nary will be broken apart by tidal forces if its center of 
mass approaches the SMBH within the tidal break-up 
radius, rbt, which is larger than the Schwarzschild radius 
of the central black hole: rgc ~ 0.04 AU. If the compact 
binary does not merge due to the perturbations by the 
SMBH, it will be disrupted before the two stellar black 
holes separately inspiral and merge with the central ob- 
ject. In this case, each of the two inspiral events cannot 
be treated as an isolated two-body problem and the mu- 
tual perturbations between the stellar black holes could 
produce a gravitational wave signal in a LISA-like de- 
tector which would differ significantly from the idealized 
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waveform obtained in the case of an isolated extreme- 
mass ratio inspiral (jAmaro-Seoane et al.ll201^ . 

It is interesting to explore the effect of changing outer 
and inner orbital eccentricities on the binary merger 
time. In the upper panel of Figure [B] each curve cor- 
responds to a fixed periapsis distance and inclination 
/ = 85°. As expected from Equation 1^, larger external 
eccentricities lead initially to faster mergers, but, as the 
orbital period of the external orbit becomes comparable 
to the timescale for relativistic precession in the inner 
binary, the merger time starts to increase rapidly with 
the orbital radius, or apoapsis, of the external orbit. For 
oi — 0.1 AU, setting Tpcr w 20 AU, the binary merger 
time is 10~^ x Tq when rapo ~ 0.01 pc. This is also 
an estimate of the limiting galactocentric radius within 
which we expect Kozai oscillations to significantly affect 
the dynamical evolution of the black hole binaries. The 
lower panel of the same figure gives Tm as a function of 
ei for / = 85°, ai = 0.05 AU and 03 = 130, 100 and 
70 AU. Note that, Tq is also a function of ei, for eccen- 
tricity larger than about ei ^0.7 the effect of the Kozai 
resonance on the merger time becomes progressively less 
important since ei approaches (or exceeds) the value of 
the maximum eccentricity that would be attained by the 
binary via Kozai processes, which is mainly fixed by the 
initial inclination (emax — \/l — 5/3cos^/, neglecting ad- 
ditional sources of apsidal precession). 

4.3. Merger fraction and eccentricity distribution 

The number of GW sources from binary mergers in 
galactic nuclei A^gw : is the number of binaries residing in 
some region close to the SMBH where they can efficiently 

merge (r < Rmerger)- 
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where Ni,{r < Rmerge) is the number of stars in the en- 
closed region, /;,„ is the binary fraction for the given 
stellar type discussed earlier, fcbm is the fraction of close 
enough binaries (a < amerge) that survive in this region 
and f merge is the fraction of binaries for which Kozai evo- 
lution lead to merger. In the following we estimate the 
rate of binary mergers as well as the properties of the co- 
alescing binaries, which will determine their GW signal 
characteristics. In order to do so we use a Monte-Carlo 
approach which accounts for all the relevant parameters 
in Equation ()29p . through random sampling of binaries 
with the appropriate properties. Binary properties are 
sampled from the large phase-space of the various bi- 
nary properties, which we describe below. Through the 
random sampling, we obtain a large ensemble of binaries 
and use the approximate method described in section 
14.21 to evolve each of the binaries and check whether 
it could coalesce before it evaporates due to encounters 
with other stars. For those binaries that coalesce we 
save the binary eccentricity as it enters the LIGO band; 
these are then used to determine the overall eccentric- 
ity distribution of the GW sources. We first determine 
the relative fractions of merging binaries arising from dif- 
ferent separations around the SMBH (see Table [T]), and 
then normalize the overall rates by the expected number 
of binaries residing in that region (which depends on the 
total number of stars multiplied by the binary fraction). 



In the following, we describe the method in details. 

We considered the idealized case of binaries with equal 
components of masses 10 Mq or 2 Mq corresponding 
to BH-BH and NS-NS binaries respectively. Given the 
uncertainty in the spatial distribution and properties of 
compact binaries in galactic nuclei the results presented 
in this section constitute a set of baselines for making 
predictions about the expected gravitational wave signal 
produced by such systems. 

We sampled the binaries external orbits from the dis- 
tribution of orbital elements: 



N{a,e^)dade^ = Noa^^^dade^, 



(30) 



which generates steady-state phase-space distributions 
for an isotropic density cusp near a dominating point 
mass potential. We first used the density model of Equa- 
tion ll]) with 7 = 0.5, 1.8 to determine the evaporation 
time scale of the binary systems and we assumed that 
they follow the same density distribution of background 
perturbers, solar mass stars in this case. In a further 
set of integrations, the evaporation time-scale was com- 
puted from the combined density profile of stars and 
BHs corresponding to a mass s egregated cusp near SgrA* 
(jHopman fc Alexand"ei1l2006al ): 



p(r) == p*(r) +/9BH 



r 
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(31) 



with pBH = 1 X lO"'Af0pc~'^. The functional form of 
the stellar density profile, p*(r), is that of Equation ([T| 
with 7 = 1.4. The power law index in equation ((5(J | 
was set to fi — 2 and (3 — \.b for BHs and NSs respec- 
tively. We imposed a lower limit for the periapsis of the 
external orbit of 02(1 — 62) > 4rbt based on the fact 
that at shorter distances the binary is unstable due to 
the strong perturbations from the central black hole and 
Equations (3)- (6) become a poor description of the dy- 
namical evolution of the triple system. The distribution 
of inner semi-major axes, ai, follows the binary distri- 
bution given in Figure 1. The eccentricity of the inner 
binaries, ei, were instead sampled from a thermal dis- 
tribution, i.e. N{< ei) ex e\. We adopted a uniform 
distribution in cos(/) and a random distribution in gi. 

An approximate estimate of the mer ger ti me was ob- 
tained with the technique described in iWenI ([2003,) and 
summarized below. Including the contribution of rela- 
tivistic precession, but in the absence of gravitational 
radiation, the maximum eccentricity in a Kozai cycle 
can be estimated by considering that the Hamiltonian 
of the system, in its quadrupole form, is a conserved 
quantity. The total Hamiltonian \s T-L = kW where k = 
3Gm nTOiTO2a?/8MiQ2(l-e^)^/^ and (jMiller fc HamiltonI 
[2OOI ) 



W{gi,ei) = -2(1 - ei) + (1 - ei)cos'l 
-|-5eisin2,gi(cos^/ — 1) 
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with the last term accounting for relativistic precession. 
The system starts from initial eccentricity ei, argument 
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of periapsis gi and mutual inclination / and evolves to 
a maximum eccentricity emax and critical gcrit and Icrit- 
These latter quantities can be related to each other set- 
ting dei/dt = in equation (|22|) . this gives: 
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Using equation ([26l) to relate /crit to Cmax, and assuming 
that Oi, 02 and 62 are constant quantities, one obtains 
Cmax solving the implicit equation: 



W{gi,ei) - W^(5crit, e,„ax) = 
An approximation to the merger time is 
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Given a certain distribution of orbital elements we then 
estimated the number of binaries that would merge in 
the galactic nucleus due to Kozai resonances induced by 
the central black hole. 

Table [T] gives the total fraction of merging binaries 
for which the external orbits are sampled inside some 
fiducial galactocentric distances. We find that for BH- 
BH binaries, the SMBH-Kozai induced mergers repre- 
sent > 40% of the total number of mergers occurring at 
r < 0.01 pc. The overall fraction of mergers in models 
characterized by a cusp in the distribution of background 
perturbers (i.e., mass-segregated and cusp models) is 
slightly lower with respect to that of models with a core 
due to the shorter evaporation time-scale at small radii 
in these former systems. When we sampled the external 
orbits within larger radii (> 0.1 pc), the merger proba- 
bility increases since the evaporation time becomes large 
with respect to the merger time of the binaries. How- 
ever, at larger radii, both evaporation processes and GR 
precession become more important in quenching Kozai 
cycles that contribute only ~ 10% to the overall merger 
rate. In the table we also give the total fraction (/') 
of binary systems with merger time in isolation shorter 
than their evaporation time-scale. These systems would 
merge even with no help from the Kozai mechanism. 
Comparing these numbers with those corresponding to 
mergers induced by the gravitational perturbations from 
the SMBH, we see that that the process described in this 
paper could contribute only a small fraction of LIGO 
sources in galactic nuclei, nevertheless, these are qualita- 
tively different than regular mergers, as they contribute 
high eccentricity GW sources as shown below. 

The total numbers of GW sources from compact bi- 
nary mergers in Milky- Way like nuclei can now be in- 
ferred by multiplying the fractions we obtained with the 
total number of binaries in such regions, i.e. the frac- 
tion of compact objects binaries in the stellar population. 
For the cusp and core models we assume the fraction 
of COs resemble their p opulation in the field. Follow- 
ing |Ho2ma^^3Ifixander 2006b we take the fraction of 
BHs and NSs to be 0.001 and 0.01 of the total stehar 
population, respectively (in the absence of mass segre- 
gation). For the mass se gregated models, we fo l low the 
distributions obtained bv lHopman fc Alexandeij ()2006bO 
(this time after mass segregation occurred); i.e. we take 



nsHir) oc r^^ and nArs(r) ex r^^-^ radial distributions, 
with the total numbers normalized such that 1800 BHs 
and 374 NSs reside inside 0.1 pc . We assume the binary 
fractions follow our discussion in section [2] (0.1 and 0.07 
for BHs and NSs, respectively). Taking these numbers 
we calculate the total numbers of merging compact bina- 
ries in a given Milky- Way like galactic nucleus (see Table 

Figure [7] shows the number distribution of eccentric- 
ities at /gw — 10 Hz (i.e., the beginning of the LIGO 
band) for merging BH-BH binaries. This plot was ob- 
tained by integrating the set of equations ([2T1) -(|24 |) and 
sampling the initial conditions as detailed above but 
considering only systems with inclinations / > 80°. 
The permitted parameter space was defined by impos- 
ing the evaporation time of the binary to be longer 
than its merger time. As additional constrain we im- 
posed the merger time of the binary in isolation to be 
shorter than the local evaporation time (i.e., TQ{ai, ei) > 

Tev/x/I — ei) since these systems would rapidly merge 
even without undergoing Kozai resonance. 

In all cases, at /gw = 10 Hz about ^ 90% of merg- 
ing systems have ei < 0.1. However, in the cusp and 
mass-segregated distributions a small fraction (about 
^ 0.1%) of merging binaries have very large eccentric- 
ities ei ~ 1 at the moment they enter the LIGO band. 
These are cases in which the gravitational radiation ef- 
fect dominates t he evolution within one Kozai cycle (see 
e.g. Figure 8 of H^[200l . About 0.5% of the merg- 
ers have eccentricities larger than 0.5 at 10 Hz. Sys- 
tems that reach such high eccentricities at the first Kozai 
maximum will be very short lived, which decreases the 
probability that they can be observed. On the other 
hand, e ccentric mergers mig ht be detectable to larger dis- 
tances (jO'Learv et al.ll2009[ '^ 1 Gpc for advanced-LIGO 
and for 10 Mq black holes, e.g.,) and greater black hole 
masses than circular mergers. Finally, we note that the 
fraction of eccentric mergers found here is a conservative 
estimate due to the lower limit that we imposed for the 
periapsis of the external orbit (i.e., 02(1 — 62) > 4rbt). 
In fact, we expect that for smaller periapsis, the binary 
could evolve through a non-secular evolution that could 
lead to very large eccentricities in only a few years (e.g.. 
Figure HI). 

4.4. GW rates from coalescing binaries in galactic 

nuclei 

In the previous section we calculated the total num- 
ber of merging compact binaries, for a given galactic 
nucleus. The total GW rated from such mergers de- 
pends on the rate at which the binary progenitors of 
such coalescing binaries are replenished. The GW rate 
form coalescing binaries in a galactic nucleus is therefore 
Tgw — /gw X Lrep, where T^-ep is the rate of binaries 
replenished into the GC timescale (if no replenishment 
occurs, we take the replenishment rate to be the total 
number of primordial binaries divided by the age of the 
Universe) . We estimate the replenishment rate assuming 
it is dominated by diffusional replenishment (see below), 
but we discuss various other mechanisms which may con- 
tribute to increase the replenishment rate, beyond our 
basic estimates. 

Diffusional replenishment: The simplest form of re- 
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Fig. 7. — Distribution of eccentricity at /gw = 10 Hz for 
BH-BH mergers driven by Kozai effects in a Milky Way-like nu- 
cleus. The external binary orbits were gen erat ed from the distri- 
bution of orbital elements given in equation I I30I I where we adopted 
/3 = 0.5 (core; continue-line histograms), 1.8 (cusp; dashed-line his- 
tograms) and 2 (relaxed cusp; dash-dotted histograms). The dis- 
tribution of background perturbers from which we computed the 
binaries evaporation times was either the broken power-law model 
of equation JTJ with 7 = /9 for the "cor e" and the "cusp" mod- 
els or the density model of equation l lijll l for the mass-segregated 
model. The distribution of semi-major axes, ai , followed that given 
in Figure 1, while the internal eccentricities ei were sampled from 
a thermal distribution (i.e., A'^(< ei) oc e^). Lower panel gives 
the details of the eccentricity distribution for ei < 0.1. For the 
cusp and the mass-segregated models most of the mergers occur at 
ei > 10-1. 



plenishment is through diffusional migration of binaries 
outside the region of influence of the SMBH in the central 
reg ion. The timc scalc for such processes (e.g. as studied 
by IHopmanI |2009() is the relaxation timescale at which 
the system approaches a steady state with no memory 
of the initial conditions. The replenishment rate assum- 
ing such a steady state is of the order of the number of 
CO binaries, given a nucleus model, divided by the re- 
laxation time, and we get Fqw = -^gw/ min (r^ep, th)- 
The calculated GW rates, given these values are shown 
in Table 1. The relaxation timescales are taken from Fig. 
[21 these are approximately constant for the cusp model 
throughout the central pc {Trep{r) ^ Const = 10^ yrs), 
and are longer than a Hubble time for the core models 
(and we therefore take Fqw = Nq-^/th)- 



Comparison with other GW sources and rate 
estimates: When compared with other GW sources 
outside galactic nuclei the GW rates from the processes 
discussed here are t ypically lower than most other es- 
timates discussed in lAbadie et al.l (J201 CD: reflecting the 
small number of compact objects residing in a galactic 
nucleus compared with the total number of such objects 
in a galaxy. The maximum BH-BH coalescence rate we 
find (for the segregated cusp model), provides a coales- 
cence ra te which is fiv e times larger than the low esti- 
mates in lAbadie et al.l |2010), 10 % of the realistic rate 
and 0.1 % of the high/optimistic rate estimate. The 
contribution to NS-NS coalescence rate is small in all 
cases; the highest rate we find (for the cusp model) is 
10 , 1 and 0.1 % of the low, realistic and high estimates 
in lAbadie et al.l (|2010D . respectively. The overall contri- 
bution of GW sources from galactic nuclei is therefore 
likely to be small, however, the eccentric signature of a 
fraction of these sources is qualitatively different than 
regular GW sources. 

Other replenishment processes: Various other 
processes may affect the replenishment rate of compact 
binaries and/or their progenitors in the GC. In the fol- 
lowing we qualitatively discuss such processes, but we 
do not account for them in our rates estimates (beside 
star formation, effectively integrated in into the assump- 
tions used in the observations based estimates), as their 
complete analysis is beyond the scope of this paper. 

Triple disruption: IPeretsI (2009; see also Ginsburg and 
Perets 2011) suggested that disruption of triple stars 
could leave behind a binary in a close orbit around the 
SMBH; this could serve as a continuous source of replen- 
ished binaries close to the SMBH. This study estimated 
a triple disruption rate which is ~ /triple = 3 — 5 % of 
the binary disruption rate of massive stars (OB stars; 
we do not consider disruption of triple compact objects 
as these are likely to be rare). In half of these cases a 
binary is captured around the SMBH, where its compan- 
ion is ejected. If the S-stars observed close to the SMBH 
in the GC arise from a binary disruptions and/or sim- 
ilarly the observed population of HVSs arise from such 
processes, then these could be used to calibrate the bi- 
nary and hen ce triple disru ption rates. C apture rate cal- 
culat i ons by iPerets et al.l (^2007: see also iMadigan et al.l 
I2009t IPerets et al.ll2009[ ). suggest that the binary dis- 
ruption scenario (induced by massive perturbers scatter- 
ing) could be consistent with the observed number of 
the young B-stars in the GC (and similar processes are 
likely to be at work in other galaxies; IPerets fc Alexandeii 
120081) . If this is the case then binaries are continuously 
replenished into close regions around the SMBH (and 
into eccentric orbits around the SMBH), and our use of 
the observed number of B-stars could be justified, with 
some important differences. Only ~ 1.5 — 3 % of the 
observed B-stars would be the result of triple disrup- 
tion and would actually be binaries. On the other hand 
the captured binaries are likely to have smaller separa- 
tions than typical binaries in the field, which are more 
likely to form compact binary progenitors that would ef- 
ficiently merge. Taken together, it is likely that such 
a replenishment mechanism would produce a GW rate 
which is 10-20 times lower then the high estimate for the 
NS stars. This scenario, however, is not likely to pro- 
duce progenitors massive enough to form SBHS, such as 
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TABLE 1 

GW RATE ESTIMATES AND THE TOTAL FRACTION OF BH-BH AND NS-NS BINARIES MERGING DUE TO KOZAI RESONANCES IN A GALACTIC 
NUCLEUS (/), AND FRACTION OF MERGERS THAT WOULD OCCUR IN A GALACTIC NUCLEUS WITH NO HELP FROM THE CENTRAL SMBH (/'). 

All THE GW rate estimates are per Milky- Way like galaxy. 



the observed young 0/WR stars in the GC (which reside 
in a disk like configuration and have low eccentricities). 
GW rate from BHs coalescence are therefore not likely 
to gain much from this process. 

In-situ star-formation: As mentioned above, our GC 
contains a large population of young massive 0-stars, 
many of which reside in a stellar disk. These stars were 
most likely to form in-situ following the fragmentation of 
a gaseous disk form ed from an in-falling gaseous clump 
(jGenzel et al.|[20Tol and references therein). Such stars 
formation burst may occur continuously throughout the 
evolution of the stellar cusp around the SMBH. As men- 
tioned earlier eclipsing and close binaries are observed 
among the GC young stars, suggesting star formation as 
an additional process which can repopulate the binary 
population of the GC. Continuous formation of binaries 
in the GC is part of the assumptions made in our high 
GW rate estimate. 

The other assumptions in the high estimate are related 
to the binary population properties and their distribu- 
tion around the SMBH. Although some pioneering work 
have bee n done on the propert ies of binaries formed in 
the GC ([Alexander et al.ll2008[) . these are generally un- 
known. Another important issue is whether the orbits of 
binaries formed in a disk could later dynamically evolve 
to achieve close approach to the SMBH, where efficient 
secular evolution as studied here can occur. Although 
the 0-stars observed in the disks are not observed closer 
than 0.05 pc to the SMBH, various mechanisms were sug- 
gested for their inward migration closer to the SMBH. 
IBaruteau et al.l (|2011[ ) discussed the disk migration of 
binaries through similar migration scenario of planets in 
gaseous disks; they show that the binaries do not only 
migrate but also shrink their inner orbit and possibly 
merge before they finish their migration. Madigan et al. 
(|2009t 120 111) proposed that binaries (and single stars) 
could be excited to high eccentricities through an ec- 
centric disk instability scenario, and their periapsis ap- 
proach could become small enough for the binaries to be 
disrupted. Both these mechanisms suggest the existence 
of binaries at close approach to the SMBH. In the disk 
migration scenario the binaries' orbit around the SMBH 
are likely to be relatively circular, and the binaries loose 
their orbital energy and migrate on a short timescale 
of the order Tmig ~ 10^ yrs; it is not clear, however, 
whether the gaseous disk enabling the migration would 
extend close to the SMBH, and/or allow for migration in 



the innermost regiontQ. In the eccentric disk instability 
binaries remain at large semi-major axis orbits around 
the SMBH, but gradually increase their eccentricity, and 
therefore come closer to the SMBH only at periapsis. 
The timescale for the close approach is relatively short, 
Tecc ~ /ew X 10^, and binaries may be left at eccentric 
orbits at the end of this process. 

Cluster infall: Nuclear stellar clusters may result 
from the continuous infall of globular clusters that dis- 
perse close or in the close vicinity of the SMBH (see 
lAntonini et al.l[2011bl and references therein). In partic- 
ular such clusters may harbor an inner core cluster of 
BHs that formed during the cluster evolution. If these 
BHs are retained in the cluster this mechanism may con- 
tribute to the BHs population in the GC. Further work 
needs to be done to explore the implications of this sce- 
nario to the replenishment of binaries in the GC. 

Resonant relaxation: In our discussion on replenish- 
ment processes above we considered the replenishment 
of binaries through non-resonant relaxation processes 
in which stars diffuse through random encounters with 
other stars. Resonant relaxation could change the eccen- 
tricities and inclination of stars much faster than regu- 
lar relaxation. One may therefore consider the timescale 
Trep for diffusional replenishment to be the scalar reso- 
nant relaxation timescale Trep = Trr, which could be 
much shorter than the non-resonant relaxation timescale 
assumed in the low estimate (see Fig. [2]) ; this could raise 
the low estimates by 1-2 orders of magnitude. For the 
high estimate, resonant relaxation may support the as- 
sumption that the distribution of stars rapidly achieve 
a steady state thermal distribution of eccentricities, as 
assumed in our models. 

The role of the much faster VRR, which changes the 
inclination of an orbit around the SMBH is less clear. 
It depends on the dynamical effect on the inner binary 
inclination in respect to its orbit around the SMBH. We 
may consider two extreme possibilities. If the binary 
orbit conserves its mutual inclination, i.e. acts as a stiff 
object, than VRR would not affect the secular evolution 
of the inner binary. However, if the binary orbit plays 
the role of a gyroscope, an initial co-planer configuration 
(e.g. zero mutual inclination), in which Kozai evolution 

^ It is interesting to note, however, that compact binary mergers 
in a gaseous disk could provide additional and different channel 
for the formation of GW sources, if compact binaries exist in the 
gaseous disk, as could be the case of long lived AGN disks. We 
leave further study of these issues to future work. 
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is ineffective, could be transformed into a high mutual 
inclination configuration at which Kozai evolution plays 
an important role. To further complicate the problem, 
the timescales for VRR could be short enough as to be 
comparable to the Kozai timescale, in which case the 
coupling of both processes should be considered. All of 
the above mentioned issues for VRR of binaries require a 
dedicated study which is beyond the scope of this paper. 

5. COMBINED KOZAI-CYCLES AND TIDAL 
FRICTION EVOLUTION OF STELLAR 
BINARIES NEAR SMBHS 

In the previous sections we discussed the secular evo- 
lution of CO binaries near a SMBH. A similar secular 
evolution could be important for stellar binaries and/or 
MS star-CO binaries. In these cases GW emission which 
served to dissipate energy and induce the inspiral of the 
CO binaries is not likely to play an important role. How- 
ever, binary stars can dissipate energy through tidal fric- 
tion. The secular evolution of the triple system, SMBH- 
stellar binary, can therefore lead to a similar outcomes 
to those stu died in the context of stellar triples (see also 
|Perets"20 09D. 

[Mazeh fc Shahainl (|1979| ) first suggested that Kozai cy- 
cles and tidal friction in triples can induce their inspi- 
ral an d the formation of compact binaries. Later stud- 
ies by lEggleton fc Kisseleva-EggletonI ()2006[ ) suggested 
that short period binaries typically form through this 
process (which could be even more effective than previ- 
ously though__whcn acco unting f or octupole approxima- 
tion; [Naozieral. 2011b). iPerets fc Fabrycky (2009) sug- 
gested that such process can induce mergers of binary 
stars and the formation of blue stragglers. Such mergers 
could be therefo re a possible source of young stars at the 
Galactic center (jPerets fc Fabrvckvll2009l : lAntonini et al.l 
l2011bt ) . Similarly other types of close binaries and their 
products could be formed in a similar way. We may 
therefore expect that the perturbative effect of SMBHs 
on binaries provide an efficient channel for the formation 
of close binaries and their products, such as X-ray bina- 
ries as well supernovae and gamma ray bursts (following 
the merger of two WDs or NS, or through the formation 



of a WD accreting binary; see also lThompson|[201 ID . We 
note, that the binary fractions of MS stars are gener- 
ally higher than those of COs, as discussed in section 
2, and tidal fraction is more efficient in inducing binary 
shrinkage/coalescence than is GW emission. Taken to- 
gether, we therefore expect the rate of shrinkage/merger 
of MS stars to higher than our estimtes of CO binary 
mergers. Since this study focuses on secular evolution of 
CO binaries, we postpone futher quantitative discussion 
of the evolution of binaries with MS star components to 
another study. 

6. CONCLUSION 

In this study we explored the secular evolution of com- 
pact binaries orbiting a massive black hole, and consid- 
ered the effects of the stellar environment of the binaries. 
We have shown that the SMBH can drive the binary (in- 
ner) orbit into high eccentricities, at which point the BHs 
can efficiently coalesce through GW emission, at times 
shorter that the evaporation times of the binaries in the 
hostile stellar environment of the SMBH. Such coalescing 
binaries can have non-negligible eccentricity when they 
enter the LIGO band, and therefore have a unique ob- 
servable signature. In addition we found that binaries at 
very close orbits to the SMBH can interact with it on 
dynamical timescale and present peculiar evolution lead- 
ing to several GW pulses over a short timescale. The 
latter sources are likely to be rare, and not be detected 
with current or planned instruments. Finally, we note 
that similar processes can be important for the forma- 
tion and merger of close stellar binaries, in which cases 
the GW dissipation leading to a the binary inspiral is 
replaced with tidal friction. 
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